Clustering instability in a freely falling granular jet 
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This paper investigates a clustering instability of a freely falling granular jet composed of 100/im 
glass spheres. The granular flow out of a circular nozzle starts out spatially uniform and then, 
further downstream, breaks up into well defined clusters. The role of air is investigated in this 
phenomenon by changing the ambient air pressure down to 1/5000*^ atm. An optical method is 
used that measures inhomogeneities in the flow in order to quantify the growth of the clusters. 
Clustering is observed down to the lowest pressure and the presence of air leads to larger drops but 
does not initiate the drop formation. The analysis shows that the drop size is set by fluctuations 
on the order of the size of the particles at the nozzle. 

PACS numbers: 45.70.-n, 83.60.Wc, 47.50.Gj 
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I. INTRODUCTION 

Granular media often appear to behave like ordinary 
fluids. One can pour sand into a bucket or use it in hour 
glasses. This paper investigates how a freely falling gran- 
ular jet emanating from an aperture becomes inhomoge- 
neous and starts to form droplets (Fig. ^ similar to an 
ordinary fluid column breaking up due to the Rayleigh- 
Plateau instability ^J] . Despite this apparent similarity in 
behavior there are considerable differences between fluid 
and granular flows. The fluid jet instability is driven by 
the surface tension of the liquid. Dry, non-cohesive gran- 
ular media, however, do not possess surface tension, so it 
is surprising to observe clustering in a granular jet. 

Due to the lack of surface tension, the clustering insta- 
bility must be driven by something else. In general, in- 
homogeneities in granular flows are quite common due to 
friction between particles and geometrical constraints im- 
posed by boundaries. This can lead to arching and jam- 
ming. In the case studied here, however, the clusters form 
in the absence of boundaries. This distinguishes it from 
other clustering phenomena in granular flows such as den- 
sity waves in funnels and vertical pipes |2) bl IJj IM Id 13 ■ 

The clustering of freely falling granular jets appears to 
be generic. Recent experiments revealed that after the 
impact of a large sphere on a loosely packed bed of small 
particles a surprisingly tall granular jet emerges [a M, 
Hol . In this case as well the jet breaks up into droplets 
whose diameter is comparable to the jet diameter. This 
is another example of drop formation in the absence of 
nearby boundaries. 

This study investigates, as a function of air pressure, 
drop formation from the flow out of a funnel. Air is im- 
portant when dealing with small grains. Viscous drag 
can easily exceed the particle's weight for the 100 ^m 
glass spheres used in this experiment. Moreover, inves- 
tigating the role of air enables a comparison with recent 
experiments on drop formation in an underwater granu- 
lar jet |lllll2l |. In my experiment I can tunc the influence 



of the surrounding "liquid" - namely the air. Studying 
the effect of air in this phenomenon is important for pin- 
pointing the clustering mechanism. In particular, I study 
the properties of the drops, such as their size and growth 
as they freely fall, as a function of pressure. 

The typical drop size decreases with decreasing pres- 
sure down to some threshold pressure below which the 
size remains unchanged. Thus, air influences the clus- 
tering but is not necessary for initiating the process. At 
all pressures the drops grow as they fall. The growth 
is well fit with a gravitational stretching function which 
indicates that perturbations of the order of a grain size 
at the nozzle serve as nucleation points for the drops. 
At higher pressures the jet atomizes below some depth. 
The depth at which disintegration occurs increases with 
decreasing pressure. At the lowest available pressure, 
p « 1/5000"^ atm, the jet does not disintegrate within 
the experimentally accessible range of depth. 

The following section gives an overview of inhomo- 
geneities that can arise in granular flows. Section III 
describes the experimental method. The next section 
contains the experimental results and is followed by a 
discussion and conclusion in sections V and VI. 



II. BACKGROUND 

Intermittency and clogging are intrinsic features of 
granular flows. This is due to friction between grains 
and geometric constraints that prevent grains from flow- 
ing past each other when the flow is confined by bound- 
aries (jamming). For finer particles air can also in- 
duce intermittency. This is observed in various sys- 
terna, such as flow down a vertical pipe or in hourglasses 
0, 0, IM 0) 0) ll^' In my experiment the system un- 
dergoes a change from a dense granular funnel flow to a 
freely falling flow. Both regimes are susceptible to inter- 
mittency and clustering. In the following, different mech- 
anisms that create inhomogeneities in granular flows and 
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FIG. 1: 100 micron glass spheres draining from the 4 mm wide circular nozzle at atmospheric pressure. This porous funnel is 
the same as used in the experiment. Pictures are taken at different depths z below the nozzle, (a) at the nozzle; (b) z = 25 
cm; (c) z = 55 cm; (d) z = 75 cm; (e) z = 130 cm. 



gases are reviewed. 



A. Density waves 



Density fluctuations can propagate through a flow, 
such as in vertical pipes |3j 14 ^^ funnel flow. Granu- 
lar funnel flow and the related hopper flow have been 
extensively studied in the literature. The flow out of a 
nozzle can become inhomogeneous in a 2D as well as a 3D 
funnel system |2, Il4| . As the grains move downward un- 
der the influence of gravity, arches form near the nozzle 
that can temporarily hold up the flow before they col- 
lapse. The presence of these density waves is enhanced 
when the grains are rough and is suppressed when the 
grains are smooth spheres |3|. 



granular flow out of an aperture can become oscillatory 
[a, i7, 13] ■ This so-called "ticking" is due to pressure gra- 
dients inside the bed that create a back-flow of air into 
the nozzle. 

In cases where the interstitial fluid is a liq uid, interest- 
ing flow instabilities have been observed [ill Il2l | . Nicolas 
|l2l | studied a system similar to mine, namely a granular 
jet emanating from a nozzle, but in the presence of a liq- 
uid. Depending on the grain and liquid parameters, the 
jet can either remain homogeneous or become unstable 
and form blobs similar to the clusters in my experiment. 
The origin of this instability still remains unclear. Nico- 
las showed that treating the suspension jet as a fluid with 
some effective viscosity cannot explain this instability. 



B. Interstitial fluid effects 



C. Inelastic Clustering 



Apart from grain-grain and grain-boundary interac- 
tions, interstitial fluid such as air can profoundly affect 
the flow. This is especially important for systems such 
as mine where the particle size is well below 1 mm [l^. 
In this regime hydrodynamic forces can influence the dy- 
namics strongly via viscous drag and pressure gradients 
building up inside the bed. The latter occurs because the 
movement of air is impeded inside a porous medium. In 
typical experimental settings (f ^ 1 m/s, p ~ 1 atm) this 
effect does not matter much for grain sizes > 1 mm, but 
becomes important for smaller media. As a result, the 



A granular gas is dissipative due to the inelastic nature 
of the collisions between the grains. Therefore, without 
external energy su pply the gas cools down with time and 
eventually freezes |la,ll3jll3- However, it does not cool 
homogeneously, but forms clusters. When a fluctuation 
increases the density locally, the collision rate goes up in 
this region. Due to the increased dissipation, the granular 
temperature decreases which in turn lowers the pressure. 
The resulting pressure gradient will enhance migration of 
particles into that region, thereby increasing the density 
even further. Eventually, these regions grow into clusters. 



D. Cohesion 

Humidity in the air and surface charges can induce co- 
hesive forces between particles [l5|. Condensates from 
ambient humidity create hquid bridges causing particles 
to stick together. As with interstitial fluid effects, the in- 
fluence of cohesion depends strongly on the size and the 
density of the particles. Smaller and lighter media are 
more susceptible to cohesive forces. Cohesion does not 
cause intermittency per se, but can dramatically change 
the rheological properties of granular media. In extreme 
cases it leads to clumping and caking. These effects can 
be controlled but not completely eliminated. In this ex- 
periment, I keep the relative humidity at 50%, which 
provides enough ions to eliminate surface charges, but 
does not result in clumping. However, small residual co- 
hesive forces could potentially induce clustering in the 
freely falling jet. 



III. EXPERIMENTAL METHOD 

In this paper I exclusively study d = 100 /xm spheri- 
cal glass beads (Mo-Sci Corp., p = 2500 g/1) (Fig. [ij. 
In order to investigate the clustering in the emerging jet 
quantitatively, an optical system was set up to measure 
inhomogeneities in the jet - similar to ones used in previ- 
ous studies for measuring intermittent granular flows in 
hourglasses and vertical pipes 0) la, 01 • The basic idea is 
to measure the light intensity of a laser shining through 
the falling grains. The fluctuations in the intensity allow 
me to record inhomogeneities in the flow. This method 
does not measure density variations, but rather the emer- 
gence of undulations on the jet surface and gaps. 

The schematics of the setup are shown in Fig. [21 The 
funnel is mounted inside and near the top of a 14 cm 
wide acrylic tube whose diameter is large enough so that 
the particles are always far from the wall. A 9 cm wide 
cylindrical reservoir feeds the grains into a porous tube. 
At the bottom of this tube a disc with a 4 mm circu- 
lar aperture is attached from which the particles emerge 
and then freely fall down the tube. There is a remote 
controlled shutter beneath the nozzle to initiate and stop 
the flow. The pressure inside the acrylic tube can be 
pumped down to 0.02 kPa. For air, this corresponds to 
a mean free path of 0.3 mm (« 3 d) [19j. The pressure 
is monitored with a pressure gauge mounted on the lid 
(Granville-Phillips, Convectron gauge 375). An oil filter 
(K.J. Lesker micromaze foreline trap) and a shut-off valve 
was installed between the system and the vacuum pump 
to avoid contamination of the system with oil vapor that 
might emanate from the vacuum pump. 

A prerequisite for this experiment is to ensure a steady 
non-changing flow out of the nozzle. It is known that air 
pressure gradients across the bed can cause oscillatory 
flow ("ticking") out of the nozzle [l3j- Even equalizing 



the pressure in the reservoir with the pressure at the noz- 
zle is not sufficient for steady flow conditions. To elimi- 
nate ticking, I made the tube coming from the reservoir 
porous. It is made out of metal mesh (0.0014" wire mesh 
325x325) that is permeable to air, but not to the parti- 
cles, thereby allowing pressure equalization through the 
boundaries (Fig. (SJ. This ensures a stable flow rate at 
all pressures (Fig. 0}. 

The whole optical system is mounted on a platform 
that is moveable in the vertical direction and can be 
lowered to 2.1 m below the nozzle. Since the width of 
the freely falling particle flow exceeds the diameter of 
the laser beam (5 mW laser diodes from z-bolt.com), the 
light is spread out to a sheet to capture the entire hor- 
izontal spread of the flow. This is achieved by shining 
the light through a 8 mm diameter glass rod. The beam 
is also focused in the horizontal plane by a lens to en- 
sure that the beam waist is smaller than the emerging 
structures. The vertical width is 0.3 mm (« 3d). After 
the laser sheet passes through the particle flow, it is re- 
focused by a cylindrical lens onto a photodiode (Silonex 
SLD-68HL1D) that measures the intensity. A current 
amp is used to convert the diode current, which is lin- 
early proportional to the incident light intensity, into a 
voltage signal. The signal is inverted since the quantity of 
interest is the amount of light that is screened by the jet. 
We refer to this signal as the blockage B. The frequency 
response is flat up to at least 100 kHz. The signal is 
recorded with an A/D converter card at a sampling rate 
of 48 kHz. The baseline of the signal changes with verti- 
cal position due to imperfections of the acrylic tube. The 
signal is scaled with the baseline at each position before 
each experimental run. 

In order to convert time scales into length scales, the 
local average velocity of the flow is measured by cross- 
correlating the intensity signal of two closely spaced laser 
sheets (vertical distance= 5.1 mm). In order to avoid 
cross-talk between the two photodiodes through scat- 
tered light a shutter is mounted near the cylindrical lens. 
A 2-channel spectrum analyzer (SRS SR780) obtains the 
cross correlation between the two signals in real time. 

When calculating autocorrelations, the mean of the 
signal is subtracted and the autocorrelation is normal- 
ized to 1 at zero time delay At. In that way, the auto- 
correlation approaches zero at large At. Near the nozzle, 
the signal to noise is low and a noise floor appears in the 
autocorrelation. The electric noise floor is measured di- 
rectly by obtaining the autocorrelation of the laser signal 
without the granular jet. It is then subtracted from the 
autocorrelation. 

Cohesion between particles is caused by liquid bridges 
and/or electrostatic interaction due to charge buildup 
|l5l |. The latter is especially significant in dry atmo- 
spheres. To ensure stable conditions the laboratory is 
controlled at 50% relative humidity. Moreover, the walls 
of the reservoir are grounded to avoid build up of charge 
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FIG. 2; Setup. The granular jet emerges from the reservoir 
at the top of the 14 cm wide acryhc tube. The tube is con- 
nected to a vacuum pump and the pressure is measured with 
a pressure gauge mounted on the Ud. The optical platform 
can be adjusted to any height between the nozzle and 2.1 m 
below. The light from a 5 mW (A = 650 nm) laser diode is 
focused by a lens, then passed through a glass rod to spread it 
out and finally refocussed with a cylindrical lens onto a pho- 
todiode. An identical laser sheet (not shown) is mounted 5.1 
mm below the first sheet. The signals of the two photodiodes 
are passed to two current amps that also invert the signals. 
The two signals are cross-correlated with a spectrum analyzer 
and one of them is sent to an A/D converter card. 



through friction. After each run, the beads are exposed 
to ionized air to neutralize any charges that might have 
built up during the run. 



IV. RESULTS 

In order to compare the clustering at different pres- 
sures we need to know the flow rate out of the nozzle and 
the velocity of the jet as a function of pressure. As men- 
tioned above, a porous tube is used to ensure steady flow 
at higher pressures. Furthermore, care has been taken to 
prepare the glass spheres in the reservoir in the same way 
after each run so that the packing fraction remains con- 
stant. Figure 01 shows the volume flow rate at different 
pressures. This was measured by recording how long it 
takes for the full reservoir to drain 158 cm^/s at a time. I 
do not let the reservoir empty completely, since the flow 
rate changes when the fill height is low. The flow rate 
is nearly constant for all pressures: At atmospheric pres- 
sures the flow rate is 1.82 cm^/s, which is slightly lower 
than at p = 0.04 kPa, where it is 1.87 cm'^/s. This few 
percent difference is negligible for the subsequent analy- 
sis. The constant flow rate allows a direct comparison of 



FIG. 3: Schematics of the reservoir. The reservoir consists 
of a 9 cm wide acrylic tube which feeds a 16 mm wide porous 
tube. At the bottom of the porous tube I attached a metal 
disc with a 4 mm circular aperture. The inside is covered 
with metal mesh and is grounded. The reservoir sits on a ring 
that has several holes in it to allow for pressure equalization 
between the reservoir and the lower part of the tube. The 
depth z is measured from the nozzle. 



the clustering at different pressures. 

In order to convert the signals recorded in the time 
domain into the length domain, we need to know the ve- 
locity at each height. Figure shows the velocity of the 
jet versus depth at four pressures. In the absence of any 
hydrodynamic drag force the velocity v should just fol- 
low V = \/2gz + Vq, where z is the distance as measured 
from the aperture, vq the velocity at the nozzle and g is 
the acceleration of gravity. When {v^ — Vq)/2 is plotted 
against the depth z, the resulting curve is linear for sim- 
ple free fall with a slope equal to the acceleration a = g. 
vo is constant within a few percent at all pressures. This 
is consistent with the constant flow rate found earlier. 
Down to the lowest available pressure we find a linear 
relationship between (w^ — Vq)/2 and z with a slope close 
to or equal to g. At higher pressures deviations are ob- 
served when the jet starts to disintegrate into a cloud of 
particles (Fig. [Tl(e)). Whenp = 101 kPa (= 1 atm), this 
happens around 1.2 m below the nozzle (Fig. |5l(d)). At 
this point the mean velocity stops growing with depth 
and the previously sharp cross-correlation peak becomes 
broad. When the pressure is lowered, the disintegration 
starts further downstream (Fig. |S1 (c)). At the lowest 
pressure (Fig. EKa)), no disintegration is observed in the 
experimentally available range of depth. 

The low drag on the jet at p = 101 kPa is surprising 
given that the viscous drag on a single grain is substantial 
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FIG. 4: Volume flow vs. time. The solid and dashed lines 
are linear fits with respective slopes of 1.87 and 1.82 cm^/s. 
(■), p = 101 kPa; (►), p = 62 kPa; (a), p = 26 kPa; (v), 
p = 3 kPa; (0), P = 0.04 kPa. 



and would lead to a terminal velocity of 0.76 m/s accord- 
ing to the Stokes formula for viscous drag on a sphere: 
-Fstokes = ^TTfJidv. A possible explanation is that the air is 
essentially trapped inside the jet so that it appears solid. 
Treating the jet as a porous medium, one can estimate 
the time it takes for air to penetrate the jet. Inside a 
porous medium, the pressure obeys a diffusion equation 
|20|. The diffusion constant \s D = {Pok)/{n{l - (/))), 
where Pq is the ambient gas pressure, /i the dynamic 
viscosity of air, and k the permeability of the granular 
medium at a packing fraction cj). The permeability is 
an empirical constant that is well approximated by the 
Carman-Kozeny relation ^: k = d'^{l - (l)f / (ISOcf)'^) . 
The value for D is 0.3 m^/s after substituting numerical 
values for the constants: fi — 1.8 • 10^^ Pa s, Pq — 101 
kPa and (j) = 0.5. The latter is a typical packing fraction 
for a random loose pack of spheres. This value decreases 
as the jet falls and gets stretched. The typical time r for 



air to diffuse into the jet is therefore t 



'jet 



/{2D) = 6.7 



/iS. This is small compared to other typical timescales of 
the system, such as the time it takes for a grain to fall 
its own diameter: d/v « 10"** s. This means that the 
jet is permeable to air. Therefore, a more sophisticated 
hydrodynamic description is needed to explain the low 
drag on the jet at atmospheric pressure. 

Figures|niand[7|display typical time traces of the signal 
at p — 0.027 kPa and atmospheric pressure, respectively. 
Just below the nozzle, z = 2 cm, drops have not yet 
formed and the signal fluctuates only slightly. At z = 50 
cm the fluctuations have visibly increased and finally, at 
z = 150 cm, the signal contains clear peaks which show 
the presence of well defined drops. At atmospheric pres- 
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FIG. 5: Velocity profile of jet at different pressures. In all 
panels, vq = 0.36 m/s. The solid lines are linear fits to the 
data. The slope represents the acceleration a. (a) p = 0.027 
kPa, a = 9.8 m/s; (b) p = 0.67fcPa, a = 9.8 m/s;(c) p = 49 
kPa, a = 9.8 m/s;(d) p = 101 kPa, a = 9.7 m/s. 



sure the behavior is similar. Comparing the time traces 
at z = 50 cm, it is apparent that the drops are bigger 
at atmospheric pressure. In atmosphere only depths up 
to ss 120 cm can be probed, since the jet disintegrates 
beyond that depth. 

In order to illustrate the evolution of the drops, the 
fluctuations of the blockage signal as a function of depth 
have been plotted (Fig. EJl. The fluctuations are just 
the standard deviation, AB, divided by the mean of the 
signal < B >. In vacuum (Fig. |Hl(a)) the fluctuations in- 
crease monotonically with depth and then saturate. Sim- 
ilarly, at atmospheric pressure the fluctuations increase, 
but then slowly decrease (Fig. |Hl(b)). 

The insets in Fig. |S1 show how < B > and AB vary 
with depth. At low pressure, the mean decreases with 
depth while the standard deviation increases. This re- 
flects the increase of undulations and gaps with increas- 
ing depth. At atmospheric pressure (Fig. |S1 (b)), the 
mean varies non-monotonically. Beyond 1 m it even ex- 
ceeds the value at small depth when the jet is still com- 
pact. The standard deviation, however, stops increasing 
far away from the nozzle. 

The reason for this marked difference in behavior is 
that in the presence of air, the jet does not stay colli- 
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FIG. 6: Blockage at different depths, z, ai p = 0.027 kPa. 
The length of the arrow equals the dip position of the corre- 
sponding autocorrelation. 
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FIG. 7: Blockage at different depths, z, at p ^ 101 kPa. The 
length of the arrow equals the dip position of the correspond- 
ing autocorrelation. 



mated far away from the nozzle as it does at low pres- 
sures. This leads to an increased blockage as the jet starts 
to spread. Furthermore, there is more spray, presumably 
caused by advection of particles by the surrounding air. 
The optical signal is sensitive to this spray since it scat- 
ters the light. This explains why the mean blockage starts 
to rise again below 60 cm or so. It should also be noted 
that even if the jet is still homogeneous, the blockage 
does not reach 100% at both pressures, since the laser 
sheet is wider than the jet diameter. 

One way to analyze the signal is to choose a threshold 
and convert the blockage signal into a binary sequence as 
done by Raafat and coworkers j5|. Anything above the 
threshold is considered a drop, anything below a gap. 
The resulting histogram for my experiment is shown in 
FigEl This was obtained atp — 0.027 kPa for two heights 
and the threshold was taken to be the signal average. 
The drop histogram displays a clear peak at z = 150 
cm, while it is less pronounced at z = 50 cm. The gap 
histogram remains flat and falls off at large gap sizes in 
both cases. The high occurrence of very small drop and 
gap sizes is due to electronic noise and particle spray. 
The histograms bear some resemblance to the ones in 
the vertical pipe flow studied by Raafat et al. |5|. Fig. 
|51 shows the emergence of a typical drop size, while the 



gap size distribution is broad. 

This analysis is prone to noise. The peak in the drop 
histogram appears only below 50 cm, even though inho- 
mogeneities appear earlier (by visual inspection with a 
stroboscope) . A more sensitive analysis is the autocorre- 
lation as shown in FigllOl 

Autocorrelations are displayed for two different pres- 
sures. At both pressures, a dip and a peak in the autocor- 
relation develops, though the peak is far less pronounced 
at p = 101 kPa. The dip position does not change with 
depth at p = 0.027 kPa. At atmospheric pressure, the 
dip remains constant as a function of depth initially, but 
then moves to higher At below z = 80 cm. Also, the dip 
position a.t p = 101 kPa is significantly larger than at 
low pressure. At low depths (z < 80 cm) it is 7.0 • 10** s 
compared to 3.8 • lO""* s at p = 0.027 kPa. 

A priori it is not clear whether the dip represents the 
time scale of drops or gaps. Looking at the time traces in 
Figs. |niand[7| drops and gaps have similar sizes. The his- 
togram (Fig. ISJ confirms this observation. The drop and 
gap size distribution are both broad and fall off around 
0.008 s. Therefore, the dip reflects the typical size of 
the fluctuations in the time domain. The length of the 
arrows in Figs. and represent the dips in the respec- 
tive autocorrelations of the time traces. They show good 
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FIG. 8: Blockage fluctuations as a function of depth, (a) 
p — 0.027 kPa, (b) p = 101 kPa. Insets: Mean and standard 
deviation of the blockage as a function of depth at the respec- 
tive pressure. The error bar delineates the standard deviation 
AB. 
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FIG. 9: Histogram of drops and gaps at p = 0.027 kPa. In 
both panels the solid line corresponds to 2: = 150 cm and 
the dashed line to z = 50 cm. The threshold to discriminate 
between drop and gaps is the average of the signal, (a) drop 
histogram; (b) gap histogram 



agreement with the typical fluctuation size by visual in- 
spection. 

Summarizing the above, we identify the dip position 
with the typical fluctuation size (in the time domain) and 
find that with increasing depth it stays constant for low 
pressures and becomes larger for atmospheric pressure. 
It follows that the clusters and gaps continually grow as 
they are stretched by the increasing velocity as a function 
of depth. 

In order to understand the growth of the drops it 
is instructive to consider the case where the jet is just 
stretched due to gravity ignoring all other interactions. 
The equations of motion for two grains, one right at the 
nozzle starting with velocity vq, the other a distance Aq 
below, are respectively: 



1 2 
z = vot+ -gt 

z' = Xo+vit+-gf 



(1) 
(2) 



the initial grain separation Aq will grow in time as X{t) 
z' — z. 

Parameterizing t with z we obtain 



A(z) = Ao + ^ 
9 



/2.9A0 
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When Ao <C TT- this reduces to 



A(.)^A„J^ + l. 



(4) 



Far away from the nozzle, when z 3> #, this can be 
further simplified to 



A(z) = Ao^ 



/2g£ 



(5) 



Since \/2gz is just the velocity v(z) at depth z for z ^ -^^ 
equation [SI can be written as 



where vi = ■\/2gAo + Vq is the velocity the grain attains 
by falling a distance Aq . Due to gravitational acceleration 
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FIG. 10: Autocorrelation at diflerent heights. Left panels 
show data at p = 0.027 kPa, the right panels at p = 101 
kPa. The depth at which the signal was recorded is indicated 
in each panel. The insets in the upper panels show all the 
autocorrelations for the respective pressure. 



where T is the time it takes to fah a distance \{z) at 

2 

depth z. This time is constant for z ^ |2.. We can now 
compare T with the dip position of the autocorrelation 
which is a measure of how long it takes for a drop to faU. 
At p = 0.027, the dip position is constant at aU depths 
to within a few percent: 3.8 • 10~* s; the clusters just 
get stretched by faUing in gravity. The velocity at the 
nozzle is uo = 0.36 m/s. Therefore, Ao = 1.4 • 10^'^m — 

2 

lAd (this also justifies the approximation Aq <C 5^ to 
simplify equation OJ. At atmospheric pressure the dip 
position is larger and grows below z = 80 cm. Using 
the value of the dip position for smaller depths, we find 
An = 2.5 • lO'-^m = 2.5d. 
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FIG. 11: Dip length vs. depth at different pressures. Each 
dip position has been converted into a length scale by multi- 
plying it with the local velocity. The solid lines are fits from 
equation |3] The error bars are the size of the symbols, (a) 
p = 0.027 kPa; (b) p = 0.67fcPa; (c) p = 49 kPa; (d) p = 101 
kPa. 



In order to see how the drop growth compares with 
pure gravitational stretching, I plot the dip length as a 
function of depth measured at four different pressures 
fFig lllll . The dip positions have been converted into 
length scales by multiplying each of them with the lo- 
cal velocity that is obtained from Fig. I plot the 
gravitational stretch equation (eqn. OJ for each data set 
using the value for Aq we found from our previous consid- 
erations, so there are no free fitting parameters. The dip 
lengths have been rescaled with the particle diameter. At 
p = 0.027 kPa, the data is well fit by this equation. Un- 
fortunately, the optical setup is not sensitive enough to 
pick up inhomogeneities above 10 cm depth. The fit for 
higher pressures agrees with the data until the onset of 
atomization. At that point, the dip position grows larger 
than it would just with gravitational stretching. 

In Fig. ^]l track the dip position in the autocorrela- 
tion as a function of pressure at constant depth z = 20 
cm below the nozzle. The dip position decreases with 
decreasing pressure until about p « 0.1 kPa at which it 
stays constant down to the lowest available pressure 0.02 
kPa. As we have seen before, the air has an apprecia- 
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FIG. 12; Dip position as a function of pressure at 2: = 20 cm. 



ble effect on the drop size. The dip changes by almost a 
factor of two. Moreover, it shows that air has no effect 
below 0.1 kPa. Below that pressure the dip position re- 
mains constant over one order of magnitude in pressure. 



V.DISCUSSION 

I have shown that the clustering of a freely falling gran- 
ular jet is influenced by the presence of air, but occurs 
even at low pressures. At the lowest available pressure, 
p = 0.02 kPa, the mean free path exceeds the particle 
size by a factor of 3, so air effects should be negligible. 
Moreover, at constant depth, the position of the first dip 
in the autocorrelation function, which is a measure of 
the typical size of the structure, stays constant below 
0.1 kPa down to 0.02 kPa (Fig^J). Therefore, air-grain 
interactions can be ruled out as the initial cause of the 
clustering. 

It is worth noting that other granular phenomena that 
depend on the ambient pressure lose their pressure depen- 
dence below « 0.1 kPa. Granular size separation in vi- 
brated beds depends strongly on pressure when the grain 
size is small {d < 500/zm) [21J, |23| . It has been found that 
below 0.1 kPa air ceases to play a role. Another granular 
phenomenon that depends on air is heaping. The sur- 
face of a bed of small particles [d < 600/im) starts to tilt 
when vertically vibrated |23| . It was found that heap- 
ing dramatically decreases below 1 kPa, since the mean 
free path of air becomes comparable to the grain size. 
This suggests that it is for this same reason that the dip 
position stops changing at low pressures in Figll2l 

This experiment does not measure density, so I cannot 
distinguish between cluster growth through agglomera- 
tion and the stretching of a cluster due to gravity. Nev- 



ertheless, the latter is always present, so it is sensible 
to compare the growth in cluster size with gravitational 
stretching. 

My results indicate that the instability arises from fluc- 
tuations on the granular level at the nozzle. At the lowest 
available vacuum, where air effects are negligible, gravi- 
tation is the only external force acting on the particles. 
Indeed, the fit for gravitational stretching is good at these 
low pressures. Extrapolating the drop growth back to the 
nozzle, I find the initial size to be of the order of a grain 
diameter. 

At higher pressures deviations from the gravitational 
stretching are observed. The extrapolation yields an ini- 
tial size that is almost a factor of two larger than at low 
pressures. Moreover, the jet starts to disintegrate into a 
cloud below some depth. This is presumably due to hy- 
drodynamic interactions, since we do not observe this at 
low pressures, at least in the observable range of depths. 

I did not observe clustering for grain sizes larger than 
200/im. Visual inspection by strobing granular jets of 
larger particles did not reveal any inhomogeneities visible 
by eye. Smaller grains, on the other hand, give rise to 
strong clustering. 

It is still unclear how these grain-sized fluctuations 
grow into drops of several particles. The fact that cluster- 
ing is not observed for glass spheres > 200/im suggests 
that cohesive forces might be responsible. However, if 
cohesive forces were dominant, so that the particles just 
stick together, the drops should not grow and remain con- 
stant in size once they are formed. On the other hand, 
the velocity correlations induced by inelastic collisions of 
equal particles do not depend on the mass of the parti- 
cles and therefore should not depend on their size. They 
just depend on the coefficient of restitution. However, 
there are indications [2J, |23 that the coefficient of resti- 
tution is size dependent. Experimentally 24] and theo- 
retically [23 it was shown that smaller spheres give rise 
to a smaller coefficient of restitution. Therefore, a granu- 
lar gas of small spheres should form clusters more rapidly 
than larger ones. Since I do not have any information on 
the coefficient of restitution of the glass spheres used in 
my experiment, it is unclear whether this size dependence 
could explain the absence of clusters in granular jets of 
large particles. 

Regardless of the clustering mechanism, an initial per- 
turbation is required that acts as a seed for the cluster. 
Since gravity cannot initiate fluctuations and air-grain in- 
teractions have been ruled out, this perturbation is likely 
to come from the nozzle. This is consistent with the 
finding that the cluster size extrapolates to a few grain 
diameters at the nozzle, which is a typical fluctuation in 
granular flows. 
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VI. CONCLUSION 

In contrast to ordinary fluids, granular flows lack sur- 
face tension and are discrete in nature. Despite these 
fundamental differences, a freely falling jet becomes un- 
stable and forms drops in both cases. However, the phys- 
ical origin of the instability in the granular case is quite 
different from its fluid counterpart. 

The length scale, after tracing it back to the nozzle 
using the gravitational stretch equation, turns out to be 
of the order of a grain size at all pressures. This suggests 
that granular fluctuations at the nozzle set the cluster 
size. Moreover, I have shown that the surrounding liquid, 
in my case air, is not required for drop formation. It does 
however, change the length scales by almost a factor of 2 
and ultimately leads to the atomization of the jet. The 
latter is not observed at low pressures. 

These results might have implications on the cluster- 
ing phenomenon found in suspensions |llLll2| . They have 
typically been ascribed to an effective surface tension be- 
tween the fluid and the sediment. My results show that 
fluctuations on the granular level can propagate down- 
stream and thereby impose a length scale. This suggests 
that the granular nature of the sediment cannot be ig- 
nored and a continuous medium description is inadequate 
to explain this phenomenon. 

This study also leaves some open questions: Cluster- 
ing is not observed for glass spheres larger than 200/^m. 
This observation should give insight into the clustering 
mechanism which still remains unclear. How do fluctua- 
tions of single grains grow into clusters of many particles? 
Further studies are needed to answer this question. 

The clustering of a freely falling granular jet is a novel 
granular instability that is initiated by fluctuations at 
the nozzle that are of the order of a grain diameter. 
These fluctuations grow into clusters downstream. The 
detailed structure of the drops may reveal information 
about these fluctuations that are experimentally difficult 
to access otherwise. 
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